Large-scale bottleneck effect in two-dimensional turbulence 
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Abstract. The bottleneck phenomenon in three-dimensional turbulence is generally associated 
with the dissipation range of the energy spectrum. In the present work, it is shown by using a two- 
point closure theory, that in two-dimensional turbulence it is possible to observe a bottleneck at the 
large scales, due to the effect of friction on the inverse energy cascade. This large-scale bottleneck 
is directly related to the process of energy condensation, the pile up of energy at wavenumbers 
corresponding to the domain size. The link between the use of friction and the creation of space- 
filling structures is discussed and it is concluded that the careless use of hypo-friction might reduce 
the inertial range of the energy spectrum. 



I. INTRODUCTION 

In high Reynolds number three-dimensional turbulence 
at large wavenumbers fc, the compensated energy spec- 
trum k 5 / 3 E(k) often shows a little positive bump, be- 
fore it drops down rapidly in the dissipation range. This 
pre-dissipative bump region was first observed in closure 
computations [l| , but no physical explanation was given. 
The bump was also observed in experimental work 0, Q 
and direct numerical simulations .4]. A physical expla- 
nation of the effect was first given in the work by Her- 
ring et al. They explain that, because the wavenum- 
bers above the Kolmogorov scale are damped by viscos- 
ity, some of the nonlocal interactions are damped, which 
causes the transfer of the energy flux to be less effective, 
leading to a pile-up of energy. A more detailed investi- 
gation of this effect was performed by Falkovich Q , who 
predicted its wavenumber-dependence, and who named 
it the bottleneck. In that work it is also inferred that 
the bottleneck should become more pronounced when 
hyperviscosity is used. This was illustrated in three- 
dimensional turbulenceQ and two-dimensional turbu- 
lence with a direct energy cascade H such as magne- 
tohydrodynamic turbulence. 

In the present work, we show that in two-dimensional 
turbulence it is possible to observe a bottleneck at the 
large scales, due to the effect of friction on the inverse 
energy cascade. The possibility of large-scale bottle- 
necks was advanced by Lohse and Mueller-Groeling 
This range is in three-dimensional turbulence generally 
affected by a forcing which supplies energy to the flow, 
complicating the observation of a possible large-scale bot- 
tleneck. 

An alternative interpretation of the bottleneck in 
three-dimensional turbulence was recently proposed by 
Frisch and coworkers [lOj. They relate the bottleneck 
to statistical equilibrium, and view it as an incomplete 
thermalization. Indeed, it was stated that if the Lapla- 
cian V 2 in the viscous term of the Navier-Stokes equa- 
tions is replaced by a hyperviscous term V 2q , with the 
dissipativity a tending to infinity, a Galerkin truncated 



system is obtained. In Galerkin-truncated 3D inviscid 
turbulence, energy piles up at the highest wavenumbers, 
leading to wavemodes in equipartition at the tail of the 
spectrum. The simultaneous observation of a fc~ 5 / 3 Kol- 
mogorov inertial range energy spectrum and this equipar- 
tition range, proportional to fc 2 , was recently observed in 
direct numerical simulations [11] (DNS) and closure cal- 
culations [HI. In between these two ranges a pseudo- 
dissipation range was observed, caused by nonlocal in- 
teractions between the thermal bath and the end of the 
inertial range. 

Since we are interested in the large scales in the present 
work, we discuss the influence of infrared wavenumbcr 
truncation, due to the finite domain-size. For three- 
dimensional turbulence, this truncation influences de- 
caying turbulence and increases the energy decay rate, 
if the size of the integral lengthscale becomes compara- 
ble to the domain-size [13| . Steady turbulence, forced 
at a wavenumber ki > fco, in which fco is the smallest 
wavenumber of the domain, called the fundamental, will 
establish an infrared region proportional to k 2 between 
ko and ki , independent of the exact choice of fco [3 • 

In two space-dimensions, the effect of infrared trunca- 
tion is very different. The absence of vortex stretching 
induces a double cascade with energy cascading to the 
large scales and enstrophy going to the small scales. In 
this case, the large scales are not directly influenced by 
forcing. The energy spectrum is not necessarily decreas- 
ing for k tending to zero, if no energy-sink is present at 
low fc. In an infinite domain, the energy will continue to 
cascade to smaller and smaller wavenumbers. In the case 
of a finite domain-size, the energy, cascading backwards, 
piles up at the fundamental. This condensation was pre- 
dicted by Kraichnan [l5j and observed numerically [16| . 
The condensation process was studied in detail by Smith 
and Yakhot [It], Oil , who showed that the statistics of the 
velocity field remain close to Gaussian until the conden- 
sate starts to form. At later times the statistics become 
non-Gaussian. We note that the condensate state is not 
the only possible state of finite-size two-dimensional tur- 
bulence in the absence of an infrared energy sink. In- 
deed Tran and Bowman 119, 20] showed that if the in- 



2 



verse energy cascade is weak, i. e. , if only a small part of 
the injected energy cascades to small wavenumbers, the 
rest being dissipated around the forcing scale, no con- 
densation takes place, but a gradual steepening of the 
spectrum occurs, when the cascade reaches the smallest 
wavenumbers, correponding to the box-size. This steep- 
ening results in a fc~ 3 range at late times. This scenario 
was also observed by Chertkov et al. [21]. In the present 
work we consider the case in which the energy cascade 
is strong and in which condensation takes place at the 
fundamental, when no energy sink is present. 



II. THE LINK BETWEEN HYPOFRICTION 
AND ENERGY CONDENSATION 

To avoid condensation in two-dimensional turbulence, 
large-scale friction is often used. Friction acts as an en- 
ergy sink. The influence of friction, which in principle 
acts throughout the spectrum, can be concentrated at 
the smallest wavenumbers by using hypofriction, analo- 
gous to hyperviscosity at the highest wavenumbers. The 
Navier-Stokes equations in this case become: 

d 

— + (u ■ V)u = -Vp + iA7 2 u + / -/i/c? Q (-V 2 )- a u; 
at 1 

V-m = 0, (1) 

in which u is the velocity vector, p the pressure, v the 
viscosity, / an isotropic forcing term localized in spectral 
space around wavenumber ki, /i the friction coefficient, kf 
the friction wavenumber, and we shall call a the friction- 
ality, in analogy with the dissipativity in hyperviscous 
flows [13 ■ 

In the presence of h ypo friction, an energy sink is 
present at low k (e.g. |22l . [23|). In this case a steep 
fall-off of the energy spectrum at low k can be created 
for large a. If the solution of the hypofrictional case 
converges to a Galerkin truncated case if a — > oo, in 
analogy with the work of Frisch et aZ.[l(|, a pile- up of 
energy can be expected around kf. In the case of linear 
friction (a = 0), no noticable bottleneck is observed (e.g. 
[lil|). The spectral slope d ln(E(k))/d ln(k), is a contin- 
uously decreasing function of k for fco < k <C fej with fco 
the fundamental and ki the wavenumber where energy is 
introduced by the forcing. 

In the present study, using a closure model, we will sys- 
tematically vary the order of the frictionality from a = 
to a = 128. We will investigate the appearance of bottle- 
necks, its influence on the inertial range scaling and the 
convergence of this bottleneck to an infrared wavenumber 
truncation. 

Some evidence for large-scale bottlenecks can be found 
in literature [H, [23[ . In these studies hypofriction was 
used with values 8 and 5 respectively. The compensated 
energy spectra do show a pile up of energy at the small 
wavenumbers. A link with the creation of coherent struc- 
tures is made but no systematic study of the influence 



of a on the large-scale bottleneck was reported in these 
studies. 

Our interpretation, in the light of the present study, 
is that the bottlenecks in these works are in fact incom- 
plete condensates, in analogy with the interpretation of 
Frisch et aZ.[lO|, that the bottleneck in 3D turbulence is 
an incomplete thermalization. The coherent structures 
that are observed are closely related to the structures, 
or condensates, that would be observed in a truncated 
solution. 

Even though numerical resources increase rapidly and 
two-dimensional turbulence is computationally accessible 
at high Reynolds numbers, a systematic DNS study of the 
influence of a would be very expensive. Furthermore the 
small wavenumber range of an energy spectrum converges 
more slowly to a statistically steady state than the small 
scales. It was shown that for this kind of studies, two- 
point closure is a convenient tool[T3, E3- For the case of 
two-dimensional turbulence, it was pointed out [25| that 
in the absence of coherent vortices, closure approaches 
compare well to direct numerical simulation. The inverse 
cascade is known to be free from large coherent structures 
and its statistics are close to Gaussian [26]. In direct 
numerical simulations, if friction is absent, the pile up 
of energy at the fundamental will eventually give rise to 
the creation of two counter-rotating vortices of the size 
of the domain. 

In the long-time limit, the statistics become highly 
non-Gaussian [l7j] so that the use of the model, which are 
based on an expansion about Gaussianity, will probably 
not be justified [271429} . The time before this happens is 
supposed to be reasonably well predicted by closure. It 
is this time-interval, in which pile-up of kinetic energy 
at small wavenumbers takes place, and which precedes 
the formation of coherent structures, on which we con- 
centrate in the present work. 

III. OBSERVATION OF LARGE-SCALE 
BOTTLENECKS 

We use the Eddy-Damped Quasi Normal Markovian 
(EDQNM) closure. It was used to study two-dimensional 
turbulence by Leith [3(| and Pouquet et aZ. [3l| - In this 
closure, the evolution equation for the energy spectrum 
is solved 

[d t + vk 2 + /i(fc//fc) 2Q ] E(k, t) = T NL (k, t) + F{k). (2) 
The nonlinear transfer Tnl is 

T NL (k,t) = - [[ O kpq XV [k 2 pE(p,t)E(q,t) 

-kp 2 E(q,t)E(k,t)]^-. (3) 
pq 

The forcing is a time-independent energy input of rate c, 
localized at k = ki, F(k) — c 5(k — ki). In equation ([3]), 
A is a band in p, g-space so that the three wave- vectors 
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k,p, q form a triangle, x, y, z are the cosines of the angles 
opposite to k, p, q in this triangle. The characteristic time 
Qkpq is defined as: 
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kpq 



1 - exp(-(r) k + r)p + rj q ) x t) 
r]k + r]p + Vg 



(4) 



in which rj is the eddy damping. In the present study we 
use the following expression for the damping coefficient 



Vk = M / s 2 E(s, t)dS + vk 2 + n{k f /k) 



2o 



(5) 



In expression ([5]), three contributions to the eddy damp- 
ing are present. The first two are the standard contri- 
butions introduced in the EDQNM closure. For A we 
use the value 0.4 [32j. The third contribution is an ad- 
hoc modification of the damping to take into account the 
presence of friction at large scale. A rigorous derivation of 
the modification induced by the friction as done for the 
rotating case [33j or using an improved eddy-damping 
[HJ is not attempted in the present work. The numerical 
set-up is similar to the one in [l2j . The wavenumbers are 
spaced cis kp — /lq r p with fco the fundamental and r > 1 
a constant which determines the number of wavenumbers 
per decade. The spatial resolution for the computations 
reported here is 40 wavenumbers per decade which cor- 
responds to a resolution of F w 12 wavenumbers per oc- 
tave. Simulations are reported for a = 0, 1, 2, 4, 8, 32, 128 
as well as a friction- free simulation truncated at fco = kf ■ 
For all runs but this last one, a steady state was ob- 
tained, for which the energy spectrum is shown in Fig- 
ure [TJ For the friction-free (truncated) steady 
state will not be obtained until the energy-pile up at 
the fundamental is balanced by the viscous dissipation 
at this wavenumber, which is very small in the present 
case as the Reynolds number is large. The fundamen- 
tal would then accumulate an amount of energy which is 
enormous compared to the cases with small frictionality. 
For a better comparison between the different cases, in 
the following a spectrum is shown for which the energy 
pile-up at the fundamental is comparable to the pile-up 
for a — 128. In Figure [5] we observe that for a > a 
bottleneck is present around kf w 3. The size of this 
bottleneck increases with a. For values of 1 < a < 8, 
the bottleneck has the form of a bump. For a > 8, the 
bottleneck changes its shape towards a peak and for the 
highest values of a used in the present work, the bottle- 
neck collapses with the energy spectrum of the frictionless 
truncated case. The fc -5 / 3 inertial range is influenced by 
this bottleneck, not only on the bottleneck itself, where 
the spectrum is steeper than fc -5 / 3 , but also in the re- 
gion with wavenumbers larger than the bottleneck, and 
extending for about one decade, the spectrum deviates 
and is slightly shallower than fc~ 5 / 3 . The friction- free 
truncated case, and the high-a friction cases, show much 
similarity (at small k) with the inviscid truncated 3D 
case (at high k) studied in [l2j]. A depletion of the com- 
pensated energy spectrum is present in the vicinity of 



the peak (for k > kf in the 2D case, and for k < k t h in 
the 3D case, where k t h is the wavenumber marking the 
beginning of the equipartition range) . The effect can, in 
both cases, be explained by nonlocal interactions, as was 
done in the three-dimensional case 12[. 

It is interesting to investigate whether this depleted 
region in the compensated energy spectrum corresponds 
to the emergence of a new scaling regime. Dimensional 
analysis based on the assumption of a constant flux of 
energy and a sweeping timescale related to the advection 
of the small eddies by the velocity associated with the 
kinetic energy contained in the condensate, one can de- 
duce a scaling proportional to fc~ 3 / 2 . It is observed in the 
inset of Figure [5] that a small wavenumber range shows a 
spectral exponent close to this value (1.53). This range 
is however too short to make conclusive statements. 

It was shown by Smith and Yakhot [l7| , that when the 
condensate becomes large enough, the statistics of the 
velocity field become strongly non-Gaussian. One could 
then ask, if the present closure approach, based on an as- 
sumption of maximum randomness [35j j , is still valid. To 
estimate the extend to which the closure gives results for 
the energy spectrum close to results of DNS, we compare 
our results with DNS results from literature [ll| . In that 
study the coexistence of a large peak at the fundamental 
and a fc~ 5 / 3 range for larger wave-numbers is observed. 
The value of the energy accumulated at the fundamental 
is one order of magnitude above its inertial range value. 
In Figure [2] the peak value is also one order of magnitude 
above its inertial range value. So at least qualitatively, 
the coexistence of a large peak and the A: -5 / 3 inertial 
range is observed in DNS as well as in closure. How- 
ever this validation remains qualitative and the present 
results could only be fully supported by a detailed and 
quantitative comparison with DNS. Such a comparison 
is beyond the scope of the present paper. 

When the bottleneck becomes large enough, a small 
kink is observed in the spectrum around k = 50. This 
comes from the discretization. Nonlocal interactions are 
not resolved if the ratio of the smallest wavenumber to 
the middle wavenumber of an interacting triad is smaller 
than[3~fl| a = 2 1 / F — 1. For elongated triads, the middle 
wavenumber and the longest wavenumber are approxi- 
mately of the same size. The largest wavenumber k di- 
rectly interacting with kf can then be shown to be k ss 50 
for F = 12. This clearly illustrates the nonlocal origin 
of the pseudo-friction range, created by the bottleneck. 
In principle, the effect of the truncation can be removed 
by refining the discretization. To allow nonlocal interac- 
tions over the full 4 decades of the simulation, the number 
of necessary gridpoints can be estimated to be roughly 
700 per octave which would make the simulations pro- 
hibitively expensive. Another solution would be to ex- 
plicitly compute the first-order contributions of the non- 
local interactions [3~fl ] by expanding the nonlinear transfer 
with respect to the parameter a, which is small for very 
nonlocal interactions. This procedure has the disadvan- 
tage not to conserve energy or to mix-up the different 
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orders of the expansion. The improvement expected to 
be minimal with respect to the presented results, these 
solutions will not be attempted here. 



IV. CONCLUSION 

The present results show that the bottleneck effect is 
not only observable at the small scales, but it can also be 
created, be it artificially, at the large scales. It shows that 
by increasing a, the energy spectrum tends to the con- 
densed solution predicted by Kraichnan This con- 
stitutes a logical extension of the ideas of Frisch et al. 
to the inverse cascade in two-dimensional turbulence. It 
also shows that pile-up of energy at small wavenumbers 
influences the inertial range scaling by nonlocal inter- 
actions. The present work thereby confirms the ideas by 
Sukoriansky et al. [36], H?} , that large-scale friction can in- 
fluence all scales of two-dimensional turbulence, since the 
pile-up of energy at the small scales might trigger the cre- 
ation of coherent structures. These coherent structures 
can induce a /c~ 3 scaling[2l|, [22[, which camouflates the 



inertial range. The use of hypofriction, intended to en- 
large the inertial range, can thus have the opposite effect. 

A detailed DNS study would be necessary to confirm 
the effects observed in the present study. Even at low or 
moderate resolution, trends could probably be observed 
and answer the question whether the predictions of the 
present closure, based on an expansion about Gaussian- 
ity, are valid for the flow here considered in which co- 
herent large scale structures are present. In the authors' 
opinion, in this type of flow, as is the case for many types 
of turbulent fields, the simultaneous use of Direct Numer- 
ical Simulations at low and moderate Reynolds numbers 
and of closures at high Reynolds numbers will probably 
constitute the best strategy to gain understanding of the 
physics of turbulence dynamics. Such a study is however 
left for future work. 
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FIG. 2: The same spectra as in Figure [T] compensated by fc 5 ^ 3 and, in the inset, compensated by fe 1,5 
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